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I. INTRODUCTION 

Proton-proton collisions at s/s = 900 GeV have been 
recorded by ALICE (A Large Ion Collider Experiment) at the 
Large Hadron Collider (LHC) at CERN fl]. Hadron colli- 
sions at these energies provide an opportunity to probe Quan- 
tum Chromodynamics (QCD) under extreme conditions. The 
distinguishing feature of QCD is the mechanism of color con- 
finement, the physics of which is not fully understood, due 
in part to its theoretical intractability J2]. The confinement 
mechanism has a physical scale on the order of the proton ra- 
dius and is especially important at low momentum. 

Bose-Einstein enhancement of identical-pion pairs at low 
relative momentum was first observed in pp collisions by 
Goldhaber, Goldhaber, Lee and Pais 50 years ago |3]. Since 
then, two-pion correlations have been successfully applied to 
assess the spatial scale of the emitting source in e + e~ 0], 
hadron-hadron and lepton-hadron @|, and heavy ion |@] col- 
lisions. Especially in the latter case, this technique, known as 
Hanbury Brown - Twiss (HBTVinterferometry and being 
a special case of femtoscopy Q^, LTO], has been developed into 
a precision tool to probe the dynamically-generated geometry 
of the emitting system. In particular, a first order phase tran- 
sition between the color-deconfined and -confined states was 
precluded by the observation of short timescales 0]. At the 
same time, femtoscopic measurement of bulk collective flow, 
manifesting itself via dynamical dependences of femtoscopic 
scales ("homogeneity lengths" ITllll2ll ). provided hints that a 
strongly self-interacting system was created in the collision. 
This was further corroborated by the positive correlation be- 
tween the HBT radius and the multiplicity of the event @]. 

In particle physics, overviews of femtoscopic measure- 
ments in hadron- and lepton-induced collisions 10, IH [HI] 
reveal systematics surprisingly similar to those mentioned 
above for heavy-ion collisions. Moreover, in the first direct 
comparison of femtoscopy in heavy-ion collisions at RHIC, 
and proton collisions in the same apparatus, a virtually identi- 
cal multiplicity- and momentum-dependence was reported in 
the two systems 11411 . 

A systematic program of femtoscopic measurements in 
pp and heavy-ion collisions at the LHC will shed consider- 
able light on the nature, the similarities, and the differences of 
their dynamics. With the present work, we begin this program. 



n. EXPERIMENT AND DATA ANALYSIS 

The data discussed in this article were collected in Decem- 
ber 2009, during the first stable-beam period of the LHC com- 
missioning. The two beams were at the LHC injection energy 
of 450 GeV and each had 2-4 bunches, one of them collid- 
ing at the ALICE intersection point. The bunch intensity was 
typically 5 x 10 9 protons, giving a luminosity of the order of 

1Q 26 

cm s and a rate for inelastic proton-proton collisions 
of a few Hz. 

Approximately 3 x 10 s minimum bias pp collision events 
were identified by signals measured in the forward scintilla- 



tors (VZERO) and the two layers of the Silicon Pixel Detec- 
tor (SPD) 11511 . The VZERO counters are placed on either 
side of the interaction region at z = 3.3 m and z = -0.9 m. 
They cover the region 2.8 < T| < 5.1 and —3.7 < T) < — 1.7 
and record both amplitude and time of signals produced by 
charged particles. The minimum-bias trigger required a hit in 
one of the VZERO counters or in one of the two SPD layers 
which cover the central pseudorapidity regions |r) | < 2 (inner) 
and |t|| < 1.4 (outer). The events were collected in coinci- 
dence with the signals from two beam pick-up counters, one 
on each side of the interaction region, indicating the presence 
of passing bunches. The trigger selection efficiency for inelas- 
tic collisions was estimated to be 95-97% lfl6ll . 

The VZERO counters were used also to discriminate 
against beam-gas and beam-halo events by requiring a strict 
matching between their timing signals (see Ref. [lj] for de- 
tails). This background was also rejected by exploiting the 
correlation between the number of clusters of pixels and the 
number of tracklets pointing to a reconstructed vertex. After 
these selections the fraction of background events remaining 
in the sample of events with at least one charged particle track 
was estimated to be below 0.1%. The trigger and run condi- 
tions are discussed in detail in Ref. fl^l . 

The 250 k events used in the analysis were required to have 
a primary vertex (collision position) within 10 cm of the cen- 
ter of the 5 m long Time Projection Chamber (TPC) lfl7ll . 
This provides almost uniform acceptance for particles within 
the pseudorapidity range |t|| < 0.8 for all events in the sam- 
ple. Within this sample, we have selected events based on the 
measured charged-particle multiplicity M. The three multi- 
plicity classes were M < 6, 7 < M < 11, and M > 12; about 
70% of all events were falling into the first multiplicity class. 
The tracks used in determining the multiplicity were the same 
as those used for correlation analysis (see below) except that 
particle identification cuts were not applied. The measured 
multiplicity was converted to the charged-particle pseudora- 
pidity density dV c h/dr| by normalizing it to the pseudorapid- 
ity acceptance and by correcting it for the reconstruction ef- 
ficiency and contamination. The correction factor was deter- 
mined from a Monte Carlo simulation with the PHOJET event 
generator iflU, fl9Tl and with the full description of the AL- 
ICE apparatus and is 0.71, 0.78, and 0.81, respectively, for the 
three multiplicity bins. The estimated systematic error is be- 
low 4%. The average charged-particle pseudorapidity density 
of the analyzed event sample is (dV c h/dr|)=3.6. An alternative 
method based on SPD tracklets (16] gave the same result. 

The ALICE Time Projection Chamber (TPC) d was used 
to record charged particle tracks as they leave ionization trails 
in the Ne-C02-N2 gas. The ionization electrons drift up to 
2.5 m to be measured on 159 pad rows; the position resolu- 
tion is better than 2 mm. Combined with a solenoidal mag- 
netic field of B=0.5 T this leads to a momentum resolution 
~ 1% for pions with pj < 1 GeV/c. The ALICE Inner Track- 
ing System (ITS) has also been used for tracking. It consists 
of six silicon layers, two innermost pixel detectors, two lay- 
ers of drift detectors, and two outer layers of strip detectors, 
which provide up to six space points for each track. The tracks 
used in this analysis were reconstructed using the information 
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from both the TPC (signals from at least 90 pad rows required) 
and the ITS. Separate studies have been done with TPC-only 
and ITS-only tracks, and were found to give results consistent 
with the combined ITS+TPC analysis. The tracks were re- 
quired to project back to the primary interaction vertex within 
0.2 cm (2.4 cm) in the transverse plane and 0.25 cm (3.2 cm) 
in longitudinal direction, if ITS+TPC (TPC only) information 
is used, thereby rejecting most secondary pions from weak 
decays. The pion tracks used in the correlation analysis had 
transverse momenta between 0. 15 GeV/c and 1 .0 GeV/c. 

ALICE provides excellent particle identification capability. 
In this analysis the particle identification was achieved by cor- 
relating the magnetic rigidity of a track with its specific ion- 
ization (dE/dx) in the TPC gas. The dE/dx of the TPC was 
calibrated using cosmic rays and its resolution was shown to 
be better than 5.5%, the design value. The contamination of 
the pion sample is negligible within the momentum range of 
0.25 GeV/c < p < 0.65 GeV/c. Below and above this range 
it is on the order of 5% and is caused by electrons and kaons, 
respectively. 



III. TWO-PION CORRELATION FUNCTIONS 

The two-particle correlation function is defined as the ratio 
C (q) = A (q) /B (q), where A (q) is the measured distribution 
of pair momentum difference q = P2 — Pi, and B (q) is a simi- 
lar distribution formed by using pairs of particles from differ- 
ent events (event mixing) [20]. The limited statistics available 
(520 k identical-pion pairs with qi nv < 0.5 GeV/c) allowed us 
to perform a detailed analysis only for the one-dimensional 
two-pion correlation function C{q lm ). The q\ m is, for iden- 
tical mass particles, equal to the modulus of the momentum 
difference |q| in the pair rest frame. 

The correlation functions were studied in bins of event 
multiplicity and of transverse momentum, defined as 
half of the vector sum of the two transverse momenta, 
kj = |pt,i+Pt,2|/2. During event mixing, all pion tracks from 
one event were paired with all pion tracks from another event. 
Every event was mixed with five other events with similar 
multiplicities; ten multiplicity bins were introduced for this 
purpose. The multiplicity binning improved the flatness of 
the correlation function at <7j nv > 1.5 GeV/c. Binning events 
according to their vertex position, on the other hand, had no 
effect on the correlation function and therefore was not used. 
Alternatively to event mixing, the denominator can be ob- 
tained by rotating one of the two tracks by 180° in azimuth. 
The correlation functions obtained using this technique are 
generally flatter at high qi nv than those from event mixing. 
The difference between the results obtained utilizing the two 
techniques was used in estimating the systematic errors. 

For the correlation structures measured here, with charac- 
teristic widths ~ 0.2 GeV/c, track splitting and track merging 
in the event reconstruction are small effects overall. Their im- 
pact on the results was carefully studied with the Monte Carlo 
simulation and turned out to be negligible. 

Another apparatus effect considered is the momentum res- 
olution. Momentum smearing for single particles has similar 



effect on the correlation structures in two-particle correlations 
i.e. it smears the correlation peak, making it appear lower and 
wider. We have studied this effect with the Monte-Carlo simu- 
lation of the ALICE detector and have found that for the width 
of the correlation peak expected here the effect is on the order 
of 1%. 

Fig-IUpresents two-pion correlation functions measured by 
ALICE in pp collisions at y/s = 900 GeV, as a function of 
event multiplicity and transverse momentum kj. The denom- 
inator of the correlation function was obtained via event mix- 
ing and normalized such that the numbers of true and mixed 
pairs with 0.4 GeV/c < q lm < 0.6 GeV/c were equal. The 
q mv range used for normalization was chosen to be outside of 
the Bose-Einstein peak but as close as possible to it. The nor- 
malized distributions of positive and negative pion pairs were 
added together before building the ratio of true and mixed 
pairs. The Bose-Einstein enhancement is manifest at low q\ m . 
A slight decrease of the correlation peak width is seen as mul- 
tiplicity grows. The kj dependence is less obvious because 




FIG. 1. Correlation functions for identical pions from pp colli- 
sions at ^fs = 900 GeV (full dots) and those obtained from a simu- 
lation using PHOJET (open circles). Positive and negative pion pairs 
were combined. The three columns represent collisions with differ- 
ent charged-particle multiplicities M\ the transverse momentum of 
pion pairs kj (GeV/c) increases from top to bottom. The lines going 
through the points represent the Gaussian fits discussed in the text. 
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the correlation baseline - the underlying two particle correla- 
tion without any Bose-Einstein enhancement - is systemati- 
cally changing its shape between the low and high transverse 
momenta. 

The correlation functions were fitted by a function account- 
ing for the Bose-Einstein enhancement and for the mutual 
Coulomb interaction between the two particles: 

C{q inv ) = ((1 -X)+XK{q im ) [l + exp(-Rf nv qi v )]) D(q inv ), 

(1) 

with X describing the correlation strength and Ri m being the 
Gaussian HBT radius J2lll . The factor K is the Coulomb func- 
tion integrated over a spherical source of the size 1 fm. It is at- 
tenuated by the same factor X as the Bose-Einstein peak. The 
factor D{q lw ) accounts for long-range correlations, like those 
arising from jets and/or from energy and momentum conser- 
vation, and plays an important role in the analysis as will be 
discussed later. 

Neglecting the Coulomb interaction K(qi nv ) = 1 the fit 
function reduces to 



C(qinv) = [l + X exp(-Rj m ,ql JV )] D(q inv ) 



(2) 



The difference between the /?j nv values obtained with and 
without the Coulomb correction is less than 0.05 fm. 

While the Gaussian fit captures the bulk scales of the cor- 
relation, at low q mv the data points lie above the fit line. This 
feature was observed previously in pion correlations from par- 
ticle collisions. An exponential fit 



C{qi„v) = [1 + X exp(-R inv q inv )} D(q inv ) 



(3) 



matches the data better. However, contrary to the Gaussian 
^inv, the Ri m parameter from Eq. (O does not have a straight- 
forward interpretation as the "radius of the source". We have 
used both functional forms and leave a detailed investigation 
of the correlation peak shape to future studies. In order to 
make the connection to established systematics at lower en- 
ergy particle and heavy-ion collisions, a careful treatment of 
the long-range correlations, visible as a slope in the baseline 
of the correlation developing with increasing transverse mo- 
mentum and represented by the factor D(q mv ) in Eqs. (DO, is 
crucial. 

In order to better understand the shape of the correlation 
baseline we have calculated correlation functions for pp col- 
lisions events generated by the model PHOJET and propagated 
through the ALICE detectors, performing an identical analy- 
sis for the simulated events as for the measured ones. The re- 
sults are shown as open circles in Fig.Q] The model does not 
contain the Bose-Einstein effect, hence the lack of the peak at 
low q mv is expected. At low kj and low multiplicity, the model 
predicts a flat correlation function. However, as kj increases, 
long-range correlations start becoming visible as a distortion 
of the correlation function baseline similar to that seen in the 
experimental data. 

The accuracy of our simulation in describing the correlation 
baseline was verified with unlike-sign pion pairs. The multi- 
plicity and kj dependence of the 7i + 7t~ functions is shown in 
Fig. |2] Correlation structures for non-identical pions include a 
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FIG. 2. One-dimensional correlation functions for Tt + Jt~ pairs from 
pp collisions at ^fs = 900 GeV. The columns and rows are defined 
as in Fig.Q] 



mutual Coulomb interaction peak, here limited to the first bin 
at lowest g; nv , and peaks coming from meson decays which 
should be correctly modeled in the event generator. Therefore, 
one can directly compare simulations with data. In Fig. [2] the 
simulated correlation functions agree reasonably well with the 
experimental data. This suggests that the same model (PHO- 
JET) can be used as a reasonable estimate also for identical 
particles to describe the correlation baseline under the Bose- 
Einstein peak. The presence of resonance peaks (like the Kg 
one at qi nv =412 MeV/c) and the fact that the simulated corre- 
lations for identical and non-identical pion pairs have different 
slopes, on the other hand, indicate that unlike-sign pion pairs 
cannot be directly used for the denominator of the identical 
pion correlations. 

The procedure employed to extract the HBT radii with 
Eq. (Q]) using the PHOJET baseline is as follows. First, the 
simulation points shown in Fig. Q] are fitted with the 2nd-order 
polynomial 



D(q im .) =a + bqi, 



(4) 



Subsequently, the experimental correlation function is fitted 
by Eq. fl}, taking the D(q mv ) from the PHOJET fit and adjust- 
ing X and Ri m . The two fits are represented in Fig.Q]by the 
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lines going through the simulation and experiment data points, 
respectively. 

In order to estimate the systematic error from the baseline 
determination we repeated the fitting procedure using a sim- 
ulation performed with the PYTHIA 12211 generator (version 
6.4.21, Perugia-0 (320) tune J3) instead of PHOJET. The 
HBT radii obtained in the two ways differ by up to 10%. In the 
following we use the average between them and we estimate 
the systematic error related to the baseline shape assumption 
to be half of the difference. 

It is interesting to see what happens with the radii if the 
slope of the baseline is neglected. Assuming a flat baseline 
D(qi m ) = a and treating a as the third fit parameter in Eq. ([TJ 
leads to Ri m values that are similar to those obtained with the 
PHOJET or PYTHIA baseline at low kj values but smaller by 
up to 30% at high transverse momenta. This is because the 
broad enhancement caused by long-range correlations will be 
attributed to Bose-Einstein correlations, giving rise to smaller 
radii (wider correlation function). The resulting apparent kj 
dependence will be discussed in Section M 

The Ri m obtained from the fit (the two highest multiplic- 
ity bins combined) is shown in Fig. [3] In order to reduce 
the statistical errors and to compare to other experiments, in 
the following sections of this article we analyze separately the 
multiplicity and the transverse momentum dependences. 
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IV. MULTIPLICITY DEPENDENCE OF THE HBT RADIUS 

The multiplicity dependence of the obtained HBT radius 
is shown in Fig. [4] and Table J] The analysis here was re- 
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FIG. 4. One-dimensional Gaussian HBT radius in pp collisions 
at */s = 900 GeV determined using pion pairs with kj = 0.1- 
0.55 GeV/c, {k[)= 0.32 GeV/c, and shown as a function of the 
charged-particle multiplicity at midrapidity (full dots). The shaded 
band represents the systematic errors (see text). For comparison, 
open symbols, red stars, and green filled boxes represent the data 
taken at the ISR JMl, RHIC Qj], and Tevatron (H, respectively. 



TABLE I. One-dimensional HBT radius in pp collisions at = 
900 GeV determined using pion pairs with kj = 0.1-0.55 GeV/c, 
(&t) = 0.32 GeV/c, as a function of the charged-particle pseudorapid- 
ity density at midrapidity. The radii were obtained using the Gaus- 
sian fit function defined by Eq. Q}. 



(cW ch /dr|) 


X 


^inv(fm) 


3.2 


0.386 ± 0.022 


0.874 ± 0.047 (stat.) toMl ( s y st ) 


7.7 


0.331 ± 0.023 


1.082 ± 0.068 (stat.) +°;°(® (syst.) 


11.2 


0.310 ±0.026 


1.184 ± 0.092 (stat.) +°;°^ (syst.) 



FIG. 3. Extracted HBT radius as a function of kj for low (black cir- 
cles) and high (red squares) multiplicity events. The error bars are 
statistical. The shaded bands represent the systematic errors related 
to the baseline shape assumption and to the fit range, added quadrat- 
ically. 



stricted to the first three transverse momentum bins kj = 0.1- 
0.55 GeV/c. The mean transverse momentum for pairs with 
<?inv < 0.2 GeV/c is (k T )= 0.32 GeV/c. The HBT radii were 
obtained by using PHOJET and PYTHIA to estimate the shape 
of the baseline, as explained in the previous section. The sys- 
tematic errors related to the baseline assumption reflect the 
difference between the two. The systematic error related to the 
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choice of the normalization and/or fit range was estimated to 
be 5%. An additional downward systematic error of 13-20% 
accounts for the difference between the event mixing and the 
rotation denominator techniques. The shaded area represents 
the three systematic errors added in quadrature. 

The charged-particle pseudorapidity density (d/V c h/dr|) of 
the lowest multiplicity bin was calculated excluding events 
with multiplicities M < 2 because these events do not con- 
tribute to the numerator of the correlation function. Including 
all events and including only events with at least one like-sign 
pair would shift the point by 0.8 to the left and to the right, 
respectively. 

An increase of the HBT radius with multiplicity is ob- 
served, consistent with the hadron-hadron collision system- 
atics above yfs ~ 50 GeV lfl3ll . While the average transverse 
momentum is similar in all four data sets, other aspects of the 
analysis, e.g. the average orientation of the momentum differ- 
ence vector, can differ so the trends, not the absolute values, 
should be compared. In heavy-ion collisions, this multiplicity 
dependence has been associated with the particle composition 
and overall volume of the final state system ll6l l26ll27ll or with 
final-state hadronic rescattering lf28ll . The relation observed 
in heavy-ion collision data |H, R ~ a + ^dAWdri) 1 / 3 , where 
a and b are constants, appears to be consistent with our data 
within our systematic errors. For high energy pp collisions, 
it has been suggested that a similar behavior could originate 
from final-state hadronic rescattering for short hadronization 
times l29ll . In an alternative scenario, the increase of the HBT 
radius with multiplicity results from the fact that the high mul- 
tiplicity pp events mostly come from hard parton scattering, 
and the hadronization length, i.e. the distance travelled by 
a parton before hadronization, is roughly proportional to the 
parton energy l30Tl . 

The fitted correlation strength A, is lower than unity, the 
value expected for the ideal Bose-Einstein case. One rea- 
son for this is the non-Gaussian shape of the peak, caused at 
least partially by pions from decays of short- (A, p) and long- 
lived resonances (co, r|, T|')- On the detector side, X can be 
reduced by the particle misidentification; this effect is how- 
ever small in our data sample. In ALICE, X decreases from 
0.37±0.03 to 0.32±0.03 between the lowest and the highest 
multiplicity, in close agreement with the E735 measurements 
at Tevatron ll25ll . A similar trend was observed by UA1 in 
pp collisions at y/s = 630 GeV/c l3lll : the fact that their X 
values were lower may have to do with the lack of the par- 
ticle identification and the resulting dilution of the correla- 
tion peak. In a final-state hadronic rescattering model lE^ll . a 
correlation strength dropping with multiplicity in high-energy 
pp collisions was attributed to the increased contribution from 
long-lived resonances in higher multiplicity events. 

An increase of the HBT radius with increasing particle 
multiplicity was recently reported by the CMS Collabora- 
tion for the same collision system and energy ll32ll . The 
authors fit the correlation peak by an exponential (Eq. (0). 
An analogous approach in our case (Table [n]i yields radii 
that are rather similar to the Gaussian ones (Table IJi once 
scaled down by y^S l32ll . In order to compare between the 
two experiments we perform a fit to an inclusive correla- 



tion (all multiplicities and k-j-'s). The exponential fit to the 
correlation functions obtained using event mixing and using 
rotation yields ^ inv = 1.61 ±0.07 (stat.)±0.05 (syst.) fm and 
^i nv = 1.31 ±0.05 (stat.)±0.22 (syst.) fm, respectively. This is 
in close agreement with the corresponding values quoted by 
CMS, 1 .72±0.06 fm and 1 .29±0.04 fm. 



TABLE II. One-dimensional HBT radius in pp collisions at yfs = 
900 GeV determined using pion pairs with ki = 0.1-0.55 GeV/c, 
{kj)= 0.32 GeV/c, as a function of the charged-particle pseudorapid- 
ity density at midrapidity. The radii were obtained using the expo- 
nential fit function defined by Eq. ifJJ. 



(dJVch/dri) 


X 


KinvVtt (fm) 




3.2 


0.704 ± 0.048 


0.809 ± 0.061 (stat.) +^ 


(syst.) 


7.7 


0.577 ± 0.054 


0.967 ± 0.095 (stat.) 


(syst.) 


11.2 


0.548 ± 0.051 


1.069 ±0.104 (stat.) toam 


(syst.) 



V. TRANSVERSE MOMENTUM DEPENDENCE OF THE 
HBT RADIUS 

One of the key features of the bulk system created in nu- 
clear collisions is its large collective flow. The fingerprint 
of this flow is a specific space-momentum correlation signa- 
ture, revealed in the transverse momentum dependence of the 
Gaussian HBT radius |6]. While quantitative comparison be- 
tween particle and heavy ion studies is complicated by exper- 
iments using different acceptances and techniques, a recent 
comparison of the HBT radii from pp and Au+Au collisions 
at RHIC indicates an almost identical pj dependence between 
these collision systems lfl4ll . Again, this raises the interesting 
question whether hadron collisions at the highest energies al- 
ready develop a bulk, collective behavior. 

The kj dependence of our measured HBT radius is shown 
in Fig. [5] The choice of the fitting method, which only weakly 
affects the multiplicity dependence of the HBT radius dis- 
cussed in the previous section, is of crucial importance for the 
transverse momentum dependence. Taking the baseline shape 
from the Monte Carlo leads to an HBT radius that is nearly 
independent of kj (filled black circles and red boxes for PHO- 
JET and PYTHIA, respectively). Assuming a flat baseline, on 
the other hand, results in a radius falling with kj (green stars). 
As discussed in the previous section, the experimental unlike- 
sign pion correlation functions are close to the predictions of 
PHOJET and PYTHIA and we consider using the average be- 
tween the two cases as baseline to be a reliable estimate for 
the HBT radii. 

The radii obtained in this fashion are summarized in Ta- 
bles |nl] and [TV] and shown in Fig. [6] where we compare them 
to RHIC and Tevatron data IU3I1 . Like for the multiplicity 
dependence, the systematic error band represents a quadratic 
sum of the error related to the baseline assumption (0-10%), 
the fit range (10%), and the denominator construction method 
(mixing/rotating, 7-17%). The lowest-fcx point is significantly 
below the RHIC and Tevatron results. It should be noted that 
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FIG. 5. One-dimensional Gaussian HBT radius in pp collisions at 
y/s = 900 GeV as a function of transverse momentum kj. Three fit- 
ting methods, differing by the choice of the baseline parametrization, 
are compared. 



the ALICE analysis was performed on a minimum bias event 
sample and the averaged charged-particle pseudorapidity den- 
sity is ((W c h/dr|)=3.6 while the Tevatron events are biased to 
high multiplicity, (dA r c t,/dr|)=14.4, similar to our highest mul- 
tiplicity bin. As visible in Fig. [3] the lowest-fcx point at the 
high multiplicity is at /?i nv ~ 1 -2 fen, approaching the Tevatron 
points. The STAR results, on the other hand, were obtained 
from events with (diV c h/dT|)= 4.3 i.e. similar to the ALICE 
case and thus a similar reasoning cannot explain the differ- 
ence. 

Two tests were performed to make sure that the low HBT 
radius value at low transverse momenta is not caused by ap- 
paratus effects. First, the analysis was repeated using only the 
ITS and thus reducing the low-momentum cut-off by about 



TABLE III. One-dimensional HBT radius in pp collisions at yfs = 
900 GeV as a function of the pair kj. The radii were obtained using 
the Gaussian fit function defined by Eq. 0. 



(k T ) (GeV/c) 


X 


^inv(fm) 




0.20 


0.35 ± 0.03 


1.00 ±0.06 (stat.) 


+0.10 
-0.20 


(syst 


0.32 


0.33 ± 0.03 


1.06 ±0.06 (stat.) 


+0.11 
-0.19 


(syst 


0.47 


0.30 ± 0.04 


0.99 ± 0.09 (stat.) 


+0.10 
-0.14 


(syst 


0.62 


0.35 ± 0.06 


0.99 ±0.11 (stat.) 


+0.10 
-0.13 


(syst 


0.81 


0.31 ± 0.06 


0.91 ±0.12 (stat.) 


+0.10 
-0.12 


(syst 



TABLE IV. One-dimensional HBT radius in pp collisions at y/s = 
900 GeV as a function of the pair kj. The radii were obtained using 
the exponential fit function defined by Eq. {3). 



(Jbr)(GeV/c) 


X 


R mw l V™ ( fm ) 


0.20 


0.63 ± 0.05 


0.94 ± 0.07 (stat.) +[J;^ (syst.) 


0.32 


0.58 ± 0.04 


0.93 ± 0.07 (stat.) (syst.) 


0.47 


0.55 ± 0.07 


0.92 ±0.10 (stat.) (syst.) 


0.62 


0.70 ±0.11 


0.98 ±0.14 (stat.) (syst.) 


0.81 


0.60 ±0.12 


0.90 ±0.16 (stat.) toil (syst.) 
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FIG. 6. One-dimensional Gaussian HBT radius in pp collisions at 
y/s = 900 GeV as a function of transverse momentum kj (full dots). 
The mean charged-particle multiplicity density was {dW c ij/dr|)=3.6. 
PHOJET simulation was used to determine the baseline of the cor- 
relations. Using PYTHIA and a flat baseline leads to systematic de- 
viations up and down, respectively; the related systematic errors as 
indicated by the shaded area. Stars and filled boxes represent the 
radii measured at RHIC fl4ll and Tevatron l25]l . respectively. 



50 MeV. This analysis yielded the same HBT radius which 
demonstrates that the energy loss is not an issue. Second, 
as seen in Fig. [3] the low-kj point is mostly driven down by 
the contribution of the low multiplicity events. Since the ver- 
tex resolution in these events is worse this might in principle 
deteriorate the momentum resolution and smear out the cor- 
relation function peak. In order to test this the analysis was 
performed without using the event vertex constraint for mo- 
mentum determination. The results, again, were unchanged. 
This, and the distinct peak in the unlike-sign pion corre- 
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lation functions in the low-multiplicity low-kj bin of Fig. [2] 
indicate that the momentum resolution is not spoiled in low 
multiplicity events. 

Even more important than the position of the first point, 
albeit related to it, is the question of the slope of the points 
in Fig. [6] Our measured HBT radius is practically indepen- 
dent of kj within the studied transverse momentum range. 
The slope crucially depends on the baseline shape assump- 
tion, as was shown in Fig. [5] The results from the experiments 
to which we are comparing in Fig. [6] were extracted using a 
flat background (although the STAR experiment also studied 
the effects of using other types of backgrounds for their data 
to account for the non-femtoscopic effects ITill ). Assuming 
that PHOJET and PYTHIA are correct such a procedure may 
lead to a misinterpretation of the low-q enhancement of the 
correlation function, that is coming from long-range corre- 
lations (most probably mini-jet like), as a Bose-Einstein en- 
hancement. As the impact of this may depend on the details 
of each experiment (certainly on the collision energy) we do 
not attempt to resolve this question quantitatively. However, 
we stress again the usefulness of non-identical pion correla- 
tion in constraining the correlation baseline. 



VI. SUMMARY 

In summary, ALICE has measured two-pion correlation 
functions in pp collisions at = 900 GeV at the LHC. Con- 
sistent with previous measurements of high-energy hadron- 
hadron and nuclear collisions, the extracted HBT radius i?j nv 
increases with event multiplicity. Less consistent is the rela- 
tion between R mv and the pion transverse momentum where 
the ALICE measured HBT radius in minimum bias events is 
practically constant within our errors and within the transverse 
momentum range studied. 
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